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Abstract 

We study asymptotically Lifshitz spacetimes and the constraints on flows between Lifshitz fixed 
points imposed by the null energy condition. In contrast with the relativistic holographic c-theorem, 
where the effective AdS radius, L, is monotonically decreasing in the IR, for Lifshitz backgrounds 
we find that both L and z may flow in either direction. We demonstrate this with several numerical 
examples in a phenomenological model with a massive gauge field coupled to a real scalar. 
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I. INTRODUCTION 



It is well known that energy conditions in general relativity play a crucial role in under- 
standing the global structure of spacetime. In particular, such conditions are used as inputs 
to the Penrose-Hawking singularity theorems as well as the Hawking area theorem of black 
hole thermodynamics. Of particular interest is the null energy condition, which states that 

T^kW > 0, (1) 

for any null vector field While this is a requirement imposed on the matter content of 
the theory, application of the Einstein equation 

Rfip — = (2) 

immediately converts this to a condition on the geometry 

R^k v > 0. (3) 

As an example of how the null energy condition is applied, consider an affinely 
parametrized null congruence specified by k^d^ = d/dX. The Raychaudhuri equation then 
gives 

df) 1 

^ = -^f 2 ~ w 1 " + - R^k", (4) 

where 9 is the expansion of the null congruence, the shear and the twist. (Here D 
is the dimension of spacetime.) So long as the congruence is twist-free and the null energy 
condition is satisfied, the right hand side is then non-positive, and we may conclude that 
d9/dX < 0. Moreover, in this case the inequality d9/dX < —9 2 /(D — 2) may be integrated 
to demonstrate that any negative expansion necessarily leads to the formation of caustics. 

Turning to AdS/CFT, the null energy condition plays a key role in the proof of the 
holographic c-theorem concerning flows between the UV and IR [IHlj- In particular, working 
in the Poincare patch and assuming <i-dimensional Lorentz invariance, the bulk metric may 
be parametrized as 

ds 2 d+1 = e 2A (-dt 2 + dx 2 d _ x ) + dr 2 , (5) 

where A is a function of the bulk radial coordinate r. For pure AdS, this function takes the 
form A = r/Lo, where Lq is the AdS radius. However, more generally, we may define an 
effective AdS radius L(r) along flows according to A'(r) = 1/L(r), or equivalently L{r) = 



1/A'(r), where primes denote derivatives with respect to r. This effective radius agrees with 
the true AdS radius at fixed points of the flow. 

The bulk metric ^ has the form of a domain wall solution, and the Ricci tensor is easily 
computed: 

iV = ~( A " + dA' 2 ) 9lIU , R rr = -d{A" + A' 2 ). (6) 

Choosing a null vector field k^d^ = e~ A d t + d r , the null energy condition in the bulk then 
translates to the Ricci condition (|3j) , which reads 

e- 2A R a + R rr = -(d- I) A" > 0, (7) 

or equivalently A" < 0. The holographic c-theorem immediately follows by taking A' = 
1/L, so that we are left with the inequality V > 0. This demonstrates that the effective 
AdS radius is monotonic increasing towards the UV, and furthermore allows us to define a 
corresponding monotonic c-function [IHI]. 

The general idea behind the holographic c-theorem is that the Weyl anomaly of the 
boundary field theory may be computed in the bulk dual through gravitational methods 
[HJ |5j. Since unitarity of the field theory is a crucial input to both the two-dimensional 
Zamoldchikov c-theorem [7] and the recently constructed four-dimensional a-theorem [5j , it is 
natural to expect a corresponding requirement on the matter content of the bulk theory. Such 
a requirement is naturally aligned with the null energy condition in the bulk, as violations 
of the null energy condition will lead to superluminal propagation and instabilities in the 
bulk [9HT2]. with corresponding violations in the holographic dual |13j . 

Questions about the nature of the bulk matter content become more pronounced when 
extending the holographic c-theorem to include higher derivative terms in the bulk theory 
(HH2I] . In particular, the null energy condition is no longer directly connected to the Ricci 
tensor according to (|3]), as the leading order Einstein equation will now pick up higher 
curvature corrections. Physically, higher derivative gravitational interactions may lead to 
non-unitary propagation of bulk gravitational modes, in which case one would hardly expect 
the dual field theory to be well behaved. Thus additional constraints on the higher derivative 
terms in the gravity sector must be imposed along with the null energy condition in the 
matter sector in order to obtain a well-behaved holographic c-function. 

While higher curvature bulk theories are a subject of much recent investigations, here we 
are instead interested in applying the null energy condition to bulk duals of non-relativistic 
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systems, and will restrict our focus to Einstein gravity in the bulk. In particular, holographic 
techniques have been developed for the study of strongly coupled critical points in condensed 
matter theories (see e.g. [22H21] for recent reviews). In a non-relativistic context, the time 
and space components of quantum critical systems do not necessarily exhibit the same 
scaling symmetry, and Lorentz invariance is broken into 

t -> \ z t, x -> Ax, (8) 

under scaling by A. Here z is the dynamical exponent, and Lorentz symmetry is broken 
when z ^ 1. Asymptotically, this Lifshitz scaling may be realized in the (d + l)-dimensional 
holographic dual by taking a bulk metric of the form 

ds 2 d+1 = -e 2zr ' L dt 2 + e 2r ' L dx 2 d _ x + dr 2 , (9) 

where L is the analog of the AdS radius, and parametrizes the bulk curvature. The Lifshitz 
scaling is then realized by shifts in the radial coordinate r. Note that <i-dimensional Lorentz 
invariance is restored when z — 1, in which case this metric reduces to that of the Poincare 
patch of pure AdS. 

In this paper, we will explore the implications of the null energy condition on flows 
between Lifshitz fixed points as well as flows between AdS and Lifshitz fixed points. In par- 
ticular, we generalize the relativistic c-theorem inequality V > to the Lifshitz case. Since 
time and space scaling are distinct, one may obtain two monotonic c-functions constraining 
the flows of L and z. However, the interpretation of these monotonic functions is not at 
all obvious; it is not clear whether they count degrees of freedom, and indeed their relation 
to Lifshitz scaling is obscure. Moreover, unlike in the relativistic case, we see that flows 
towards increasing and decreasing L are both allowed, so long as the dynamical exponent z 
is also allowed to flow [25J . 

In order to get a better understanding of the null energy condition and its relation to 
Lifshitz flows, we examine a phenomenological model where the Lifshitz geometry arises 
from a massive vector field coupled to a real scalar. By adjusting the scalar potential as well 
as the scalar coupling to the vector, we construct flows between different values of L as well 
as z. For concreteness, we take a simple cubic potential (so the scalar can flow between two 
critical points), although the general features of the flows do not depend on the details of 
the potential. 
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This paper is organized as follows. In section [ITJ we introduce flow functions L(r) and z(r) 
that reduce to constants L and Zq at Lifshitz fixed points. We then study the restrictions 
imposed on these functions due to the null energy condition. In section |III| we emphasize 



the distinction between Lifshitz to Lifshitz flows versus Lifshitz black holes. In section IV 
we study flows in the massive vector model coupled to a real scalar and construct several 
numeral examples with L and z flowing in various directions. Finally, in section V we 
make a connection between the present investigation of Lifshitz flows and related results 
for backgrounds that are only conformal to Lifshitz (namely those exhibiting hyperscaling 
violation). 

II. CONSTRAINTS ON LIFSHITZ FLOWS FROM THE NULL ENERGY CON- 
DITION 

At a fixed point, a Lifshitz scaling background may be written in the form (|9]), where L 
and z determine the fixed point. To study the flow of L and z, we may consider asymptotic 
Lifshitz spacetimes which arise from a generic bulk action 

S =^j d d+l X^~g{R + Matter), (10) 



where £ m atter will be taken to satisfy the null energy condition. Extending (|9j) into the 
interior of the bulk spacetime, we make the ansatz 

ds 2 d+l = -e^d^ + e^dxl^ + dr 2 . (11) 

In order to match with ^ in the UV, the functions A and B must satisfy 

A — > Zuyr/Luy and B — )■ r/Luv as r — Y oo, (12) 

where Luv and ^uv are constant values in the UV. 

When considering flows between fixed points, we would like to extend L and z away from 
their constant fixed point values. A natural way to do this is to rewrite the asymptotic 



condition (12) as 

A'(r) — > zuv/Luv and B'{r) — > 1/Luv as r — > oo. (13) 
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In particular, the derivatives A' and B' approach constants at the UV fixed point. This 
suggests that we define the L(r) and z(r) flow functions: 



L(r) = l/B\r), 

z(r) = A'{r)/B'{r). (14) 

It is clear that L(r) and z(r) agree with the constant Lifshitz radius L and critical exponent 
z at Lifshitz fixed points of the flow. 

We are now able to examine the consequences of the null energy condition. Since we 
consider a conventional Einstein-Hilbert action in the bulk, the null energy condition is 
equivalent to the Ricci condition (|3]). By taking null vectors k^d^ = e~ A d t + d r and k^d^ = 
e~ A dt + e~ B d x , we obtain two conditions [181 EE] 



R' r - R\ = (d- l)[-B" - B' 2 + A'B'] > 0, 

R x x - R\ = A" - B" + A' 2 - (d - 1)B' 2 + (d- 2) A'B' > 0. (15) 
Note that these conditions may be written as 

[e A+{d - 1)B (A' -B')]' > 0. (16) 

As a result, the two functions 

d(r) = -e {B ~ A) B', C 2 (r) = e ^(d-i)B^ _ ^ 

are non-decreasing along Lifshitz flows to the UV. The function C\ was previously obtained 
in [26] , where it was related to the monotonicity of the speed of light in the bulk, and hence 
to causality in Lifshitz holography. 

While these two c-functions appear to be the natural monotonic quantities along Lifshitz 
flows, unfortunately their relationship to L and z are not entirely obvious. Ideally, one would 



use (14) to rewrite C\ and C 2 in terms of L and z. However, the exponential factors preclude 
any straightforward interpretation. Note that, in the relativistic case where A = B, these 
functions reduce to C\ = —B' = —1/L and C 2 = 0. The C 2 function is trivial because of 
Lorentz invariance between x and t, while monotonicity of C\ implies the well known result 
that L increases monotonically towards the UV. 
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In the Lifshitz case, the functions C\ and C2 do not approach constant values at fixed 
points. Instead it is easy to see that 

d ~ -L e d-zo)r/L 0} C2 „ Z _0_zA e (zo + d-l)r/L _ (lg) 

In particular, both functions scale exponentially with r. Of course, we could directly rewrite 



the conditions (15) in terms of the effective L and z functions (14). The result is 



L'+0-1)>0, Lz'- (z- 1)2/+ (z- l)(z + d- 1) > 0. (19) 

Again, taking the relativistic case (z = 1) yields 2/ > for the first inequality. However, in 
general, the first inequality may be rewritten as 

L'>-{z-l). (20) 

As a result, this no longer leads to a restriction on the sign of V whenever the effective 
critical exponent is greater than one. Furthermore, by combining the two inequalities, we 
see that 

z' >-(z-l)(2z + d-2)/L (provided z > 1), (21) 
so in addition the flow of z can be in either direction as well. (Although this inequality is 



weaker than the individual inequalities in (19), it is nevertheless straightforward to verify 



that (19) does not preclude flows of L and z in either direction.) Finally, although the null 



energy condition does not appear restrictive for Lifshitz flows, it does yield the standard 



requirement that z > 1 at Lifshitz fixed points simply by setting 2/ = z' = in ( 19). 



III. LIFSHITZ VERSUS BLACK HOLE FLOWS 



While the null energy condition yields the two constraints (19) on Lifshitz flows, the 
physical implications of these constraints is somewhat obscure. To develop a better under- 
standing of the constraints, in the next section we will examine flows in a simple model with 
a real scalar coupled to a massive gauge field. However, before we do so, it is worth making 
a distinction between AdS black holes and Lifshitz flows. 

Since both holographic Lifshitz backgrounds and planar AdS black hole metrics can be 



written in the form ( 11 ), where the symmetry between time and space is broken, solutions to 



the bulk equations of motion by themselves cannot distinguish between the two cases. Thus 



the difference necessarily arises by imposing conditions on the solution. For flows between 
two Lifshitz fixed points (or between AdS and a Lifshitz fixed point), the bulk metric must 
flow asymptotically between two regions with well-defined scaling, while for black holes, the 
IR flow will reach a horizon and then a singularity. 



A. Schwarzschild-AdS black holes 



As an example of a black hole flow, consider for simplicity the pure Schwarzschild-AdS 
black hole, conventionally written as 

dr 21 



J* 2 - T 2 
ab d+i — -^o 



-r z fdt z + r 2 dx\_ x + 



r 2 f 



(22) 



where / = 1 — {tq/t) . Although this metric is not in the form (|9j), a simple coordinate 
transformation 



r — > r ( cosh 



2L /d 



2/d 



brings it to the form 



ds 2 d+1 = L 2 r 2 ( cosh 



2L /d 



4/d 



tanh 



2L /d 



dt + dx d _ x 1 + dr , 



(23) 



(24) 



where the horizon is at r = and the boundary is at r = oo. Using the definitions (14), we 
read off the effective quantities 



L(r) = L coth 



z(r) — 1 H — csch z 



2L /d' ~ VJ ^ ' 2^" 2L /d' 
By construction, both L and z start at an AdS fixed point in the UV 

(L,^) r -^(L ,1). 



(25) 



(26) 



However, their effective values both diverge as r approaches the black hole horizon. Of 
course, the Schwarzschild-AdS solution represents a thermal background, and not a flow to 
a non-relativistic IR fixed point. Thus this divergence is not surprising. 

It is interesting to note that, while the relativistic c-theorem implies that L is mono- 
tonically decreasing along flows to the IR, the Schwarzschild-AdS solution instead has an 
increasing flow to the IR, as shown in Fig. [TJ Since there is no matter, this flow saturates the 



inequalities n20v and (21) implied by the null energy condition. In fact, this first inequality, 



L' > — (z — 1), suggests that whenever z > 1 (which is the usual Lifshitz case), the flow 



z 

5r 




FIG. 1: The flow of (L, z) for a pure Schwarzschild-AdS4 background. The arrow points towards 
the IR (i.e. decreasing values of r), and the dotted line highlights the AdS value z = 1. 

of L tends to increase towards the IR. In order for L to decrease, the bulk matter must 
contribute enough energy density to overcome the pull from the critical exponent. This 
balance between tendency towards non-extremality versus added matter must be resolved 
in order to flow to a stable IR fixed point. When this is done, as we will see below, flows 
can then occur in directions of both increasing and decreasing L. 



B. Asymptotically Lifshitz black holes 

For actual examples of black holes in an asymptotically Lifshitz spacetime, we may con- 
sider the model of pTJ, or equivalently [28]. The latter formulation is given in terms of a 
massive vector field coupled to gravity 

S = i / dd+lx ^3 (R-V- - \m 2 A^) . (27) 

Here V and m are constants parametrizing the model. This system admits Lifshitz back- 
grounds of the form ([9]), where V and m are related to (L, z) by [251 ESI EH] 

VL 2 = -[z 2 + (d - 2)z + (d - l) 2 }, m 2 L 2 = (d - l)z, (28) 
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and with the vector field 



A = e zr ' L \j 2 ^ - l) dt. (29) 



Black hole flows of this system were considered in [29], while flows between AdS and/or 
Lifshitz fixed points were considered in [251 12Z] • The reason this simple model admits flows 



between Lifshitz fixed points is that the conditions (28) are invariant under the transforma- 
tion 

z^ (d ~ 1)2 , L->(d-l)-. (30) 

z z 

Hence, for an appropriate set of parameters V and m [corresponding to 1 < z < (d— l) 2 ], the 
system admits two Lifshitz solutions along with vacuum AdS. As shown in [22], a relevant 
deformation from the z < d — 1 Lifshitz branch can cause a flow to either the z > d — 1 
Lifshitz branch or AdS in the IR. 

As an example of the different possible holographic flows in this model, consider a four- 
dimensional bulk with VL\ = —64/9 and m 2 L^ = 8/3. These quantities are chosen so that 



the system admits two Lifshitz fixed points according to (28) in addition to a vacuum AdS 
fixed point: 



(L, z) = (^27/32^0, 1), (L, z) = (L , 4/3), (L, z) = (3L /2, 3). (31) 

Starting at z = 4/3 in the UV, this background can flow to either z = 1 or z = 3 in the IR 
[25] . The z = 1 fixed point is of course an AdS bulk, which in turn can be deformed into 
Schwarzschild-AdS. Similarly, the z = 3 Lifshitz background admits a black hole [2H] with 
a regular horizon, and with the massive vector field vanishing at the horizon. These flows 
are shown in Fig. [2} 

What this example demonstrates is that flows between fixed points can proceed in either 



direction of increasing and decreasing L and z. For z ~ 1, the two inequalities (20) and (21 ) 
favor flows where both L and z are decreasing towards the IR, as in the flow from Lifshitz 
to AdS. However, for z > 1, the right hand sides of the inequalities are negative, and there 
is plenty of room for flows in the opposite direction, as evidenced by the flow to z = 3 in 
the IR. 

In this simple example of Lifshitz backgrounds supported by a massive vector field, the 



fixed point criteria (28) restricts the critical values of (L, z) to lie on the curve 

' L 2 . (32) 



m 



d- 
10 




0.8 



1.0 



1.2 



1.4 



1.6 



1.8 



2.0 



L/Lq 



FIG. 2: The flow of (L, z) for the massive vector model with four bulk dimensions and with fixed 
points at z = 1 (AdS), z = 4/3 and z = 3. The solid line flows connect fixed points, while the 
dashed lines correspond to z = 1 (Schwarzschild-AdS) and z = 3 black hole solutions. 

As a result, flows connecting two Lifshitz fixed points in this model always proceed from 
smaller (L, z) in the UV to larger (L, z) in the IR. In particular, both L and z necessarily 
flow in the same direction. In order to examine more general Lifshitz flows, we will need to 
decouple the fixed point values of L and z. One way to do this is to introduce a variable 
effective mass for the vector field by introducing a scalar coupling. This is what we will do 
in the next section. 

IV. LIFSHITZ FLOWS IN A PHENOMENOLOGICAL MODEL 

As we have seen above, flows between Lifshitz fixed points can have a very rich structure, 
and do not appear to be strongly constrained by the null energy condition. Flows can also 
occur between AdS and Lifshitz fixed points, and between different AdS solutions, as shown 
in [25] in the context of six-dimensional M = (1, 1) gauged supergravity and its consistent 
truncation. The advantage of working in the context of supergravity is that it naturally 
provides a stringy context for the bulk dual. However, the details of the full supergravity 
may at times obscure the key features of the flows. Thus we turn to a phenomenological 
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model that nevertheless allows for a wide range of possibilities. 

In order to examine holographic flows between unrelated fixed points, we extend the 
massive vector model of [27J EE] by adding a real scalar 0, so that the bulk action takes the 
form 

S = / d d+1 xV^ (R - \d^(j> - \F^F^ - V{(j>) - W{4>)A^) . (33) 

For the moment, we leave the potential V((f>) and scalar coupling W(4>) arbitrary, although 
we assume approaches a constant, 0(r) = O , at fixed points of the flow. At such points, 



this system effectively reduces to (27) with V = V(0o) and m = 2W(4>o)- The fixed point 
values of (L, z) can then be extracted by inverting (28). 
We proceed by making a domain wall ansatz 

ds 2 d+1 = -e 2A ^dt 2 + e 2B{r) dx\_ x + dr\ 

A = e G{r) dt, <f) = 0(r). (34) 



The action (33) gives rise to four dynamical equations 



= (d - 1)[B" + B'(B' - A')} + \<P' 2 + e- 2A+2G W{<f)), 
= A" + (d — 2)B" + A'(A' - B') + ^' 2 - ^- 2A+2G G' 2 , 
= G" + G' [G' -A' + (d— 1)B'] - 2W((j)), 

= 0" + 0' [A' + {d- 1)B'] - dV/dcj) + e~ 2A+2G dW/d(j), (35) 

along with one constraint 

(d - l)B'[2A' + (d- 2)B'] + \ e - 2A+2G G' 2 - e- 2A + 2G W{$) - \$ 2 + K(0) = 0. (36) 

The exponential factors e~ 2A+2G are associated with contractions of the vector potential, 
and may be removed by defining 

H{r) = e-^W+GW. (37) 

Note that H(r) is simply the time component of the vector field in the natural vielbein basis. 
In order to focus on holographic RG flows, we may replace the metric functions A and 



B by the effective L and z functions given by (14). With this substitution, we obtain the 
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equations for (L, z, 0, H) 



L' = -(z-l) + 



2{d- 1) 



[z - l)(2z + d - 2) L 

L~ ~ + 2 



+ 2W(0)# 

2-1 



/2 



z + d-1 „ ^ 

'0+^7- 



dW 



H" 



L 

z + d-1 



H' 



d- 1 
(d 



H' + tH)\ Z \ d - 2 2W( ( j ) )H 2 
L / a — 1 



Dz 



H + 2W((p)H. 



(38) 



In these variables, the constraint equation becomes 



V(4>) - w(4>)H 2 - / 2 



2y H ' + l H 



+ 



(d- l)(2z + d- 2) 
T 2 



0. 



(39) 



At this point, it is worth pointing out that the first terms on the right hand sides of the 



L' and z' equations saturate the null energy condition inequalities (20) and (21), as they 
originate from the gravitational sector. The additional terms on the right hand side arise 
from the matter sector, and will contribute positively provided z > 1 [as already noted in 



(21)] and W(<j>) > 0. This latter condition on W{(j>) is a direct consequence of the null 
energy condition. 



A. Lifshitz fixed points 

Before proceeding to flows, we examine fixed points of this model where (L, z, <p, H) are 
all constant 

(L,z,<t>,H) -> (L ,z ,(P ,H ). (40) 



Substituting these constant values into the equations of motion (38) and the constraint 



equation (39) gives two classes of solutions. The first has vanishing vector potential, Hq = 0, 
and hence yields an AdS solution with z Q = 1 and the AdS radius L related to the critical 
value of the potential in the usual manner 

z = 1, V(</>o) = -d(d-l)/Ll, d 4 ,V{<f> ) = 0, #0 = 0. (41) 
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The second class is of Lifshitz form with 



V(cf> Q )Ll = -\zl + {d- 2)z + (d- l) 2 ] 

2W(<p )L 2 = (d - l)z , 
T 2 _ 2(^o - 1) 



z 

2(z - 1) 



(42) 



Note that the first two expressions map directly onto (28), while the value of Hq corresponds 



to (29). The final expression in (42) simply ensures that the scalar is at a critical point of 



its effective potential (which includes its coupling to the background vector). 

Flows between fixed points may be generated by perturbing by a relevant deformation. 
In order to examine what deformations are allowed, we perform a linear stability analysis 
around a given Lifshitz fixed point specified by (L , z , (f> , H ). To do so, we first rewrite 



the equations of motion (38) in first order form by splitting the second order equations for 



and H into pairs of equations for (0, </>') and (H, H'), respectively. We then perturb around 
the Lifshitz fixed point by taking 



L = L + eL, 

z = z + ez, 



p = 0o + e<; 



H = H + eH, 
H' = eH'. 



(43) 



For e < 1, the linearized equations of motion reduce to V' = AiV, where 



V = {L,z,0, cf>>,H,H>} 



(44) 



and 



M 



( ^ -1 g.H* Wl 

2(z i) -l)(z +d-2) z +d-l (zo+d-2)L jj2-\y n 
Lq d—1 1 

1 



rr 




2{z +d-2) u 
F3 -"0 










V 2 - H$W 2 



zp+d—l 
Lo 



(z -l)(d-l) U a 



z H 

2zp+d~2 TJ 


-2iWx 




\ 

z H 



1 



2(z -l)(zQ+d-2) 3zp+d-3 , 
I? L / 



(45) 

Here we have expanded the potential V((f>) and scalar coupling W(<f>) about the critical value 
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v(<t>) = v + - (f> ) + |y 2 (0 - O ) 2 + ■ ■ • , 

W{(j>) = W + Wt{(t> - O ) + \W 2 {(j> - O ) 2 + • • • . (46) 
At this point, it is worth summarizing the various parameters that enter into the matrix 



M.. The critical point determined by (42) may be thought of as having particular values of 
the Lifshitz parameters (Lq, zq). Then, at this fixed point, the constant values of V(<j)o) and 
W((j>o) become 

_ z l + (d-2)z Q + {d-lf _ (d-l)z 

VO- f2 ' W - Wf2 ■ l 4 'J 

1j q AIjq 



These quantities, however, do not directly enter into Ai. The linear terms in (46) are 



constrained by the last condition in (42) 



V l = HlW 1 , where H 2 = 2( *° 1} . (48) 

z o 

Since H is determined in terms of Zq, it cannot considered a free parameter. As a result, 
the linearized behavior at a Lifshitz fixed point depends only on the values of (L ,z ), the 
linear potential parameter W\ (or equivalently V\) and the quadratic potential parameters 
V 2 and W 2 . 

The solution to the first order equation V = MV is given by 

V(r) = ^V i e A * r , (49) 

i 

where {Aj} are the eigenvalues of Ai, with corresponding eigenvectors {V,}. Since we assume 
the asymptotic form of the metric (J9]), where r — > oo corresponds to the UV, relevant 
deformations (which induce a flow to the IR) will correspond to negative eigenvalues, namely 
< 0. If this flow terminates at a stable IR fixed point, then it will necessarily approach 
the IR fixed point along a direction or set of directions with positive eigenvalues, \f R ^ > 0. 

The eigenvalues of M. are easily determined, although their general expressions are not 
particularly illuminating. Much of the complication arises when the linear term is non- 
vanishing {i.e. when W\ ^ 0), as in this case the scalar and vector fluctuations no longer 
decouple at fixed points. We therefore restrict the presentation of explicit results to the case 
when W\ = 0. In this case, there is a single marginal mode, Ai = 0, with eigenvector 

Vx = \ H Lo,2z H ,0,0,^,o \ . (50) 
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This marginal direction corresponds to holding fixed, and moving along the curve given 



by 2WqLq = (d — l)zo in (42). Note, however, that while this deformation is marginal in 
the set of first order equations, movement along this direction will not satisfy the constraint 



(39). Thus this mode takes us outside of a given model, and in particular shifts the vacuum 
energy Vq (which is a constant of integration of the set of first order equations). 

There are three additional modes that keep <p fixed, but involve a combination of the 
metric and vector field. The first one of this set is always relevant, with A2 = — (z + d— 1)/L 
and eigenvector 

v 2 .(^_,£^±i Ho ,o,o,-l,^±^). (51) 

\z + a - 2 z + d - 2 z z L J 

The other two are paired up, with 

A 3 ,4 = - ^~ 1 (1±A), (52) 

and 

( 2-Az Q + {z Q + d-l){l^&) -4^o + (z + rf-l)(l T A) tj 

Vs ' 4 " ^ 2(z + d-2) H ° L °> 2 H ° J 

0,0,-1, (ab + d - 1)(1=FA) ), (53 ) 

where 




A^/ 1+ j VI ( r^ +1) < 54) 

Since we take A to be positive, A3 [corresponding to the top sign in (52)] will always be 
negative, corresponding to a relevant deformation. However, the behavior of A4 depends on 
the value of A. For 1 < zq < d — 1, we find A < 1, so that A4 < 0. On the other hand, we 
obtain an irrelevant deformation, A4 > 0, whenever zq > d — 1. 

The remaining two modes only involve the scalar field and its derivative (in the lin- 
earized analysis), and have eigenvalues 

A5 ' 6 = - Z0 ^o~ 1 (1±S) ' ( 55 ) 

along with eigenvectors 

V 5>6 = (0, 0, {Z ° + d ~^ {lTJ:) , V 2 - H*W 3 , 0, 0) , (56) 



16 



where 



AL 2 (V 2 - H%W, 



^TU-w- (57) 



Since we are assuming the absence of linear terms in V(4>) and W(<f>) at the critical point, 
it is clear that the combination m§ cfI = V 2 — HqW2 is simply the effective mass of 0. Note 
that, in the AdS case when z = 1, we find the expected relation between scalar mass and 
conformal dimension 

A 5 , 6 = -j- (d/2 ± yftd/W + mg iCff L§) (for z = 1). (58) 

The only modification for general zq is then the replacement d — > zq + d — 1. The behavior 
of depends on the sign of rn^ eS . For — [(zq + d— 1)/2L ] 2 < Tn^ eS < (where the lower 
bound is taken as a generalization of the Breitenlohner-Freedman bound [301 [31]), both A5 
and Ag are negative, so the scalar deformations are both relevant. However, for ml eS > 0, 
one of the deformations becomes irrelevant, corresponding to Xq > 0. 

In summary, after linearizing about a Lifshitz fixed point, we find a set of six pertur- 
bations, of which one is marginal. (This remains the case even when W\ 7^ 0.) However, 



this marginal direction is to be disregarded, as it is eliminated by the constraint (39). Of 
the remaining five modes, three are always relevant, while the other two will depend on 
the parameters of the fixed point. In the absence of a linear coupling, the modes become 
irrelevant in the following cases: 

A 4 > 0, when z > d — 1, (59) 
A 6 > 0, when m\ eS = V 2 - H$W 2 > 0. (60) 

Since "^2,3,4 have vanishing components along and </>', we label the flows induced by these 
modes as "vector field driven" 1 . On the other hand, since Vs 5 6 initiates flows along and 
$ ', we will label such flows as "scalar field driven". 



B. AdS fixed points 



Although our main interest is to examine Lifshitz fixed points, the simple model (33) 



also admits AdS fixed points given by (41). The stability analysis of such fixed points is 



1 These modes are actually combinations of vector and metric perturbations. 
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similar to that of Lifshitz fixed points. However, we cannot simply set zo = 1 in the previous 



results, as the second condition in (42) is no longer applicable when the massive vector is 
turned off. Nevertheless, the deformations can be classified in essentially the same manner. 
Again, there is a single marginal mode with Ai = and 

Vi = (1,0,0,0,0,0). (61) 

This marginal direction corresponds to shifting the cosmological constant, and is again 



removed by the constraint (39). The three modes found above involving the metric and 



vector field now split into a relevant metric deformation with A2 = —d/Lo and 

V 2 = (L ,rf,0,0,0,0), (62) 
and a pair of vector field deformations with 



d ± y/(d - 2)2 + 8W Ll 

As ' 4 " 2L (63) 

and 

v,^(o,o,o,»,- 1 ,««g). (64) 

The metric deformation generates a flow to a Schwarzschild-AdS black hole, while the two 
vector field modes have the expected conformal dimensions for a massive vector in AdS 
with an effective mass of rn\ eS = 2Wq. (Recall that the null energy condition demands 
W(<p) > 0, so that m\ eS must be non-negative.) For < m\ cS < (d — 1)/Lq, both vector 
modes are relevant, while for m\ eS > (d — 1)/Lq, the deformation corresponding to A4 
becomes irrelevant. 

Finally, the remaining two modes are those of the scalar field, with 

A 5 , 6 = -j- [d/2 ± yf(d/2)* + V 2 L 2 ^ , (65) 

and 



in agreement with (|58|), where mQ eS = V2 is just the scalar mass read off from the potential. 



As a result, we find that the deformations away from an AdS fixed point are all relevant, 
except for the two cases 

A 4 > 0, when m\ cS = 2W > (d - 1)/L 2 , (67) 
A 6 > 0, when m 2 cS = V 2 > 0, (68) 
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where the corresponding deformations become irrelevant. As above, flows to or from AdS 
fixed points may be thought of as either vector field driven or scalar field driven. Of course, 
since the vector field vanishes at AdS fixed points, a scalar field deformation may initiate 
a flow between AdS fixed points, but would not otherwise generate a non-vanishing vector 
field background. Such scalar field driven flows will thus maintain z = 1 throughout the 
flow, unless accompanied by a metric or vector perturbation. 



C. Flows between fixed points 

With the linear stability analysis out of the way, we now consider flows between fixed 
points. Such flows may be generated by starting at a UV fixed point (L\jy, zuv)> an d turning 
on a relevant deformation. This model admits five deformations away from any fixed point, 
of which three are always relevant. Hence it is always possible to construct a flow away 
from the UV. However, not all such flows will terminate in an IR fixed point. We find that 
generic deformations will flow to a singularity with both the scalar and vector fields running 
away to infinity. Furthermore, as discussed previously, flows to a black hole horizon are also 
possible. 

Assuming, however, that the flow reaches an IR fixed point with (Lir, £ir), it will ap- 
proach this fixed point in an irrelevant direction. For a Lifshitz fixed point (assuming van- 
ishing linear coupling at this point), the criteria for an irrelevant deformation are given 



by either (59), which requires Zir > d — 1, or (60), which requires mQ IR > 0. The former 



corresponds to a vector field driven flow, while the latter corresponds to a scalar field driven 



flow. For an AdS fixed point in the IR, the first criterion is replaced by (67), which requires 
™i,m> {d-l)/L\ n . 

At this point, it is worth noting that the flows of [27] in the massive vector model are 
naturally vector field driven, as this model lacks the extra scalar. In this case, Lifshitz fixed 
points in the IR must lie in the zir > d — 1 branch, as can be seen in the example of Fig. [2} 
Similarly, it can be seen that AdS fixed points in the IR can be reached for flows starting 
with 1 < z vy < (d — l) 2 . 

In the absence of a specific model, we forego a general analysis of flows between fixed 
points. However, it is instructive to provide a few examples of scalar field driven Lifshitz 
to Lifshitz flows. In particular, we may engineer such flows by choosing appropriate forms 



19 



of the scalar potential V{4>) and vector field coupling W ((/>). A natural, and fairly minimal 
choice would be to construct the potential V to have two critical points, corresponding to the 
UV and IR. We then design W so that the effective vector mass, ml eff , takes on appropriate 



values at the critical points, so as to yield the desired Lifshitz scaling according to (47). 



Although it is possible to choose a linear function for W, for simplicity we demand that 
both V and W have critical points at the same values of <fi. This allows us to avoid any 



linear shift in the potential implied by (48) 



The simplest model with two critical points involves taking cubic potentials 

V(<f>) = V + V 1 <j ) + V 2( j ) 2 + V 3( j ) 3 , 
W((f>) = W + W X (j) + W 2 (j) 2 + W 3 (f) 3 . 



(69) 



Although a cubic W((f>) will always become negative in some domain, we will always restrict 
so that W((f>) > to ensure compatibility with the null energy condition. This is a 
legitimate restriction since we are only interested in classical solutions to the equations of 
motion. The eight constants of the potentials are fixed by demanding that they admit an 
IR critical point at (f> = with (L IR , z IR ) and a UV critical point at (f> = (f> with (L uv , £uv)- 
We find, in particular 

zf R +(d-2)z m +(d-l) 2 



V 
Vi 

200 



T 2 
^IR 



0. 



'z 2 R + (d-2)z m + (d-lf 
T 2 

z? K + (d- 2)z m +(d- if 

T 2 
^IR 



+ (d - 2)^ uv + (d - 1) 



T 2 

-^uv 



'UV 



+ (d-2)z uv + (d- if 



T 2 



(70) 



and 



Wo 
W x 



(d - l)z m 



2 -^IR 



3(d 



2^0 



w 3( pl 



-^IR 



^UV 
T 2 



(d-1) 



gig. 

T 2 
-^IR 



•gUV 
T 2 



(71) 



It is worth noting, however, that this model does not necessarily admit a stable flow from 
arbitrary values of (Luvi^uv) t° {L\r,i z tr)i as relevant deformations from the UV could 
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FIG. 3: A scalar field driven flow with decreasing z and L. We have taken d = 3 and 



1, along 



with the fixed point parameters given in ( 72 ) . 



often flow to a black hole or singular geometry. In addition, there is a two-fold degeneracy 



of fixed points (both at the UV and the IR), corresponding to the map given by (30). Hence 
flows could originate or terminate in the dual fixed points from the ones that were originally 
desired. 



In order to go beyond the linearized analysis, we solve the equations of motion (38) 



numerically, integrating out from the IR to the UV along an irrelevant direction. The 
reason we integrate out from the IR is that, as indicated above, there are at most only 



two irrelevant deformations in the IR, given by the conditions (59) and (60). Since we are 



interested in scalar field driven flows (where flows between the UV and IR critical points 
of the potential), we naturally choose the direction for the deformation. An example of 
a scalar field driven flow towards decreasing z is shown in Fig. [3j where we have chosen a 
four-dimensional bulk and fixed point parameters 



(Luv^uv) = (3L /2,9/5), (L m ,z m ) = (L ,6/5). 



(72) 



Note that z lies in the range 1 < z < d — 1 for this class of flows. We see that there is some 
slight ringing as <fi flows away from 0o (where V(<f>o) is a local maximum) in the UV. 

It is also easy to obtain flows towards increasing z. One such example is shown in Fig. [4], 
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FIG. 4: A scalar field driven flow with increasing z and L. The fixed point parameters are given 



in ( 73 ) , along with d = 3 and 



1. 



where the fixed point values are given by 



(Luv, zuv) = (L , 6/5), (Lir, z m ) = (3L /2, 3). 



(73) 



The general feature of this flow is similar to the vector field driven flow from z = 4/3 to 
z = 3 shown in Fig. [2j However, we note that the effective value of L actually overshoots 
its IR fixed point value before finally turning around and reaching it at the end of the flow. 
This is a clear indication that the function L(r) is not constrained to be monotonic by the 
null energy condition. It is not obvious whether this has any physical significance, as the 
effective L and z functions are not necessarily observable along the flow. It is also possible 
to obtain flows towards increasing z, but with decreasing L. An interesting example is given 
in Fig. [5j corresponding to 



(Luv^uv) = (3L /2,6/5), (Lir, zir) = (A), 9/5). 



(74) 



Here we see that the effective values of both L and z flow away from the UV for some 
distance before turning around and reaching the IR fixed point. The ringing near the UV 
fixed point of the potential is also more pronounced, as both the scalar and vector fields 
participate in the flow. 
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FIG. 5: A scalar field driven flow with increasing z and decreasing L. We have taken d = 3 and 



1, along with the parameters given in (74) 



V. DISCUSSION 



In the relativistic case (z — 1), the consequence of the null energy condition directly leads 
to the holographic c-theorem, V > 0. However, once relativistic invariance is no longer 
required, the implications of the null energy condition are very much relaxed. Although it is 



possible to define two monotonic functions, C\ and C*2, as in (17), neither one of them serves 



the purpose of a useful c-function, as they do not approach constant values at fixed points 
(except when z — 1, in which case C\ coincides with the usual holographic c-function). 

For a Lifshitz fixed point, the more useful quantities (L, z) do not exhibit any obvious 
monotonicity. Using a cubic potential model, we have constructed flows between two Lifshitz 
fixed points that satisfy the null energy condition, and where L and z are simultaneously 
increasing or decreasing or where L is decreasing while z is increasing. So far, however, 
we have been unable to find any examples of flows with simultaneously decreasing z and 
increasing L. It is possible that this is a feature of the toy model, as the null energy 



condition in itself does not preclude such flows. Nevertheless, the inequalities (17) obviously 



place some restraints on the allowed flows, and it would be interesting to see whether there 
is a deeper reason why we have not found any flows where z decreases while L increases. 
Although we have focused on flows between scale invariant Lifshitz fixed points, there 
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has been recent interest in systems involving hyperscaling violation (parametrized by 9) in 
addition to the dynamical critical exponent z. Gravitational duals to such models have 
been realized in Einstein-Maxwell-dilaton theories (32H31], and aspects of the null energy 
condition have been investigated in [35H37] . Metrics exhibiting hyperscaling violation may 
be written in the conformally Lifshitz form 2 



ab d+l 



r 



20/(d-l) 



dt 2 d^d-i + dr 2 



,,2z 



(75) 



In order to make contact with (11), we transform r — > r^ d 1 ^ 61 (and rescale the coordinates) 
so that 

ds l +i = _ r 2(l-z(d-l)/e) dt 2 + r 2{l-{d-l)/9) dS 2_ i + ^2 (7g) 

Therefore, in a hyperscaling region, we make the identification 



A(r) 



1 



9 



log r, 



1 



9 



logr. 



(77) 



This may be contrasted with the behavior A ~ zr/L and B ~ r/L in the case where 
9 = 0. (Note, however, that this identification is not well-behaved in the limit 9 — > 0, as the 
scale-invariant metric functions behave as exponentials and not power-laws.) 

It is interesting to see what the flow functions L(r) and z(r) defined in (Jl4)) look like for 



the metric (76). We find 
L(r) = 



d- 1-9 



z{r) 



\){d 



d-l-9 



(78) 



where we should keep in mind that z(r) is an effective function that was designed to match 
the critical exponent only at pure Lifshitz fixed points, while z is the true critical exponent 



given in (76). By construction, z{r) and z coincide in the absence of hyperscaling violation 



[i.e. when 9 = 0). However, the effective function L(r) runs linearly in r, and hence cannot 
be assigned a fixed value in a hyperscaling region of any flow. 



For L(r) and z(r) given in (78), the null energy condition (19) gives rise to the inequalities 



(d - 1 - 9){{z - l){d - 1) - 9) > 0, (z - l)(z + d - 1 - 9) > 0, 



(79) 



2 We have continued to use 'relativistic AdS/CFT' notation, where the bulk is (d + l)-dimensional. Most 
of the recent hyperscaling violation literature uses d to denote the spatial dimensions of the field theory, 
so that the bulk is (d + 2)-dimensional. 
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as noted in [3 TJ 3 . It would be interesting to see whether these conditions can be extended 
along a complete flow between regions with different scaling behavior. Of course, the two 



functions C\(r) and C^fV) in (17) are valid c-functions regardless of the geometry. However, 
the usefulness of the functions depend on being able to express them in terms of 'physical' 
quantities such as 9 and z at fixed points of the flow. We have, unfortunately, not been able 



to find a suitable set of 9(r) and z(r) functions that extend the behavior of (77) beyond the 
fixed points. The main obstacle in doing so is to properly reproduce the power-law behavior 
implicit in the logr running of A(r) and B(r). 

A possible attempt at defining a c-function for flows with hyperscaling violation is to 
capture the power-law behavior e A ~ r a by removing the explicit r dependence in (77). 



This may be done by forming ratios, such as (A') 2 /A", (B 1 ) 2 / B" or A'/B'. For example, we 
may define 



B" + A'B' 



z[r) 



6{r) = {d-iy 



B" 



(80) 



B" + (B') 2 ' v ' v ' B" + (5') 2 ' 
so that z(r) and 9(r) approach the constant values z and 9$ at fixed points of the flow. 
This extension of z and 9 away from fixed points is not unique, but has the advantage that 



the first inequality in (16) yields simply 



(d - 1 - 9{r)){{z{r) - \){d - 1) - 9{r)) > 0, 



(81) 



which is identical in form to the first equation in (79), but now must hold everywhere along 



the flow. Unfortunately, the second inequality in (16) does not have a simple expression in 



terms of z(r) and 9{r) because A" in (16) has no natural counterpart in the definition (80) 



Note, also, that the inequality (81) does not take the conventional form of a gradient of a 



c-function, and hence does not suggest any direction for the flow of z or 9. Nevertheless, 
this discussion suggests that additional information may be captured from the null energy 



condition beyond just the fixed point inequalities (79). 



Finally, once again ignoring hyperscaling violation, we note that it may be possible to 
resolve the tidal singularity at the Lifshitz horizon by flowing into AdS2 x IR d_1 in the deep IR 
EST. Since the AdS 2 x 



geometry is obtained by taking B (r) = in (11), the definition 



pcZ-l 



of the flow functions L(r) and z(r) in (14) break down in this case. In particular, the flow 



Note that we must take d d+ 1 to match the notation of [37] . 
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to AdS2 x M. d ~ l is represented by 

L(r) oo, z(r) —> oo, with z(r)/L(r) = L 2 held fixed, (82) 

where L 2 is the radius of the emergent AdS2- Note that this is distinct from the Schwarzschild 
black hole flows shown in Figs. [l]and[2j as the Schwarzschild flows have z(r) ~ L 2 (r) — > oo 
as the flow approaches the black hole horizon. 

Regardless of the nature of the model, as long as it is a two- derivative theory of gravity, 
the power of the null energy condition is that it directly translates to a condition on the 
geometry, ([3]), and hence a condition on its geodesies and causal structure. Although there 
has been progress in understanding the implications of the null energy condition for higher 
derivative gravity, once additional terms enter the left-hand side of the Einstein equation, 
the null energy condition by itself no longer completely determines the Ricci tensor. As a 
result, additional conditions may be required in the gravitational sector in order to have a 
well-behaved holographic dual. It would be worthwhile to extend our results for holographic 
Lifshitz flows to the case of higher derivative gravity in the bulk, and to investigate what 
form these additional conditions may take. 
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